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ABSTRACT 

It has been proposed by Atiyah and Manton [1] that the dynamics of Skyrmions in IR'^"''^ 
may be approximated by motion on a finite dimensional manifold obtained from the moduli 
space of SU(2) Yang-Mills instantons in IR^. Motivated by this work we describe how similar 
results exist for other soliton and instanton systems. We describe in detail two examples for 
the approximation of the infinite dimensional dynamics of sine-Gordon solitons in H^"*"^ by finite 
dimensional dynamics on a manifold obtained from instanton moduli. In the first example we use 
the moduli space of (EIP^ instantons in and in the second example we use the moduli space 
of SU(2) Yang-Mills instantons in R^. The metric and potential functions on these manifolds 
are constructed and the resulting dynamics is compared with the explicit exact soliton solutions 
of the sine-Gordon theory. 
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1. Introduction 



There are numerous examples in theoretical physics of theories which possess topological 
solitons, which are extended structure solutions with a particle-like interpretation. In (1+1)- 
dimensions there are models which axe integrable at both the classical and quantum level, thereby 
allowing an explicit and exact study of dynamical soliton interactions. However, in {d + 1)- 
dimensions there are no known integrable models for c? > 1 if one restricts to lagrangian theories 
with the physically desirable property of Lorentz covariancc. We are therefore led to the study 
of non- integrable nonlinear field and gauge theories, which even at the classical level is a diffi- 
cult task. Invariably analytical and/or numerical approximations are required. For d > 1 an 
investigation of classical soliton dynamics through a numerical simulation of the field equations 
requires the use of a large amount of computing power for even a modest study. Furthermore, 
it is unclear how to address the quantum aspects of these theories from a consideration of the 
full field theory; which are in fact often non-renormalizable. It is therefore evident that there 
is sufficient motivation for attempting to find analytical approximation schemes. This paper is 
concerned with the construction, implementation and accuracy of such schemes for the situation 
in which there are forces between static solitons. 

For theories of Bogomolny type there are no forces between static solitons and the second 
order nonlinear partial differential equation which determines the static soliton solutions can be 
reduced to a first order equation. For such theories there is an analytical approximation scheme 
which works remarkably well; namely the geodesic approximation [2] . Here one assumes that the 
n-soliton configuration at any fixed time may be well approximated by a static n-soliton solution. 
The only time dependence allowed is therefore in the dynamics of the finite dimensional n-soliton 
moduli space A4n (effectively the parameter space of static n-soliton solutions). A lagrangian on 
A4n is inherited from the field theory lagrangian, but since all elements of A4n have the same 
potential energy the kinetic part of the action completely determines the dynamics. It determines 
a metric on A4n and the dynamics is given by geodesic motion on Mn with respect to this metric. 
The geodesic approximation has been applied to the study of soliton dynamics in many theories of 
Bogomolny type including BPS monopoles [3], vortices at critical coupling [4] and cr-model lumps 
[5] . For some of these systems numerical simulations of the field theory equations of motion have 
been performed and demonstrate the validity of the geodesic approximation. Quantum soliton 
interactions can also be studied by quantization of the finite dimensional system, which has been 
studied in great detail for example in the case of BPS monopoles [6]. 

Turning now to the situation in which forces exist between static solitons (ie non-Bogomolny 
theories) it is clear that if a truncation to a finite dimensional manifold is to be performed then 
a substitute for the moduli space must be found. In the special case that the theory may be 
considered as a perturbation of a Bogomolny theory then one may use the moduli space of this 
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associated theory. The difference now is that there will be both a metric and a potential function 
on A4n so that the dynamics will no longer be given by geodesic motion on Mn but by a more 
complicated dynamics determined from the lagrangian on A^„ . Such a scheme was applied to 
a planar Skyrme-like model [7], considered as a perturbation of the 0(3) cr-model (a Bogomolny 
theory), and found to produce accurate results when compared with numerical simulations of 
the field equations. Other theories which are natural perturbations of Bogomolny theories are, 
of course, non-EPS monopolcs and vortices which are not at critical coupling. Very recently this 
scheme has been applied to the study of vortex scattering at near critical coupling [8]. 

In the general case it is more difficult to find a substitute for the moduli space. One approach 
is to use the parameter space of an n-soliton configuration obtained by patching together n copies 
of a single soliton. This will certainly be an adequate manifold to describe solitons which arc 
well separated when inter-soliton forces are weak. However, it may not provide an adequate 
description of solitons which are close together, when the solitons become greatly distorted and 
lose their individual identities. In particular this patching method is inadequate for the Skyrme 
model, where the minimal energy 2-soliton can not be described by a patched configuration of 
two single solitons. One proposal for the Skyrme model is to use the union of gradient flow 
curves descending from a suitable family of static saddle point solutions [2]. Recently Atiyah 
and Manton [1] have shown that an apparently very similar manifold can be obtained from the 
moduli space of SU(2) Yang-Mills instantons in by calciilating holonomies. They have shown 
that it provides a good truncation of the 1-soliton sector and have investigated the differential 
and topological structure in the 2-soliton sector. This connection between instantons and solitons 
is still somewhat mysterious, and a clear explanation of the underlying reason for its existence 
is still lacking. There are some technical difficulties which together with the complexity of the 
calculations are responsible for the lack of a result at present on the metric and potential functions 
on this manifold. 

In this paper we describe in detail two simpler examples of an Atiyah-Manton-like procedure 
as applied to a truncation of the sine-Gordon model in (l+l)-dimensions. We show how to 
obtain suitable finite dimensional manifolds from both (CIP^ instanton moduli in and SU(2) 
Yang-Mills instanton moduli in ]R'*. By calculating the metric and potential functions on these 
manifolds we obtain an approximation to the dynamics of solitons in the one and two soliton 
sectors. These results are compared with the exact sine-Gordon soliton solutions. We find that 
the method works well; however, an unexpected discovery is that a constraint must be placed 
on the orientation of the instantons in order to avoid a spurious bound state in the two soliton 
sector. 

The plan of the paper is as follows. In section 2 we briefly review the Skyrme model and 
describe how we view the sine-Gordon theory as its lower dimensional analogue. In section 3 we 
compute sine-Gordon kink flelds from (DP^ instantons and discuss the finite dimensional manifold 



3 



this produces. Section 4 describes an alternative possibility to the previous section where we use 
Yang- Mills instantons. Section 5 is devoted to an analysis of the dynamics generated on the finite 
dimensional manifolds obtained in the previous two sections. Finally in section 6 we make some 
concluding remarks and observations. 

2. Topological Solitons in the Skyrme and Sine-Gordon theories 

The Skyrme model is a nonlinear field theory which provides an effective description of low 
energy hadron physics. Originally introduced by Skyrme [9] this model has been the subject of 
a great deal of recent study following the observations of Witten [10] . The field of the model is 
an SU(2) matrix U with lagrangian density given by 

C = -lTr(a^C/a''t/-i)+^Tr([C/-ia^C/, U-^d^UWU-^d^'U, U-^d''U])+^Tx{U -2) (2.1) 

where a;'' = (f, x), fi = 0,1,2,3 are the spacetime coordinates in Minkowski space with metric 
T]'^" = diag(— 1, 1, 1, 1) and Tr denotes trace. We have scaled away some phenomenological 
constants by a suitable choice of energy and length units, and m is a small constant equal 
to the pion mass. For the purposes of our discussion we can neglect the pion mass term and will 
therefore consider the Skyrme model with m = unless otherwise stated. 

In order to consider configurations with finite energy the boundary condition 

C/(x) ^1 as |x| ^ oo (2.2) 

is imposed. This effectively compactifies space from to .S*"^, so that at a fixed time, [/ is a 
map J7(x) : ^SU(2). Since is the group manifold of SU(2) the relevant identity is the 
homotopy group relation 

n3(5') = Z (2.3) 

which implies that to each configuration there may be associated an integer B, which is conserved 
and represents the winding number of the field as a map from space to the target manifold. It 
is known as the topological charge and counts the number of solitons in the configuration. It is 
physically identified as the baryon number so that a soliton is interpreted as a baryon. Explicitly 
B is given by 

B = ^ j eijkTr{{U-'diU){U-^djU){U-^dkU))d''^ (2.4) 

where indices range over the values 1,2,3. There is a lower bound on the energy of a given 
configuration in terms of the soliton number ie E > 12n'^\B\. 
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The static one-soliton solution (ie -B = 1) is known as the Skyrmion and has the hedgehog 

form 

[/ = exp(z/(r)x.) where x = — , r = |x| (2.5) 

and denote the Pauh matrices, /(r) is a profile function, with boundary conditions /(O) = tt , 
/(oo) = 0, which has to be determined numerically [11]. This Skyrmion is situated at the origin 

but may be moved to any position by a translation. The orientation may also be changed by 
conjugation of J7 by a fixed clement of SU(2). The Skyrmion has energy Ei = 1.23 x 127r^ and 
so does not saturate the topological lower bound on the energy. This allows the possibility of a 
bound state of two solitons in the B = 2 sector, and indeed a bound state exists, although it is 
not of the hedgehog form. The complicated nature of the 2-soliton sector and the forces between 
individual solitons make the study of soliton dynamics in the Skyrme model a difficult problem. 
This is the reason that the approximation scheme proposed in [1] is both useful and interesting. 

We now describe the sine-Gordon model in a way that illustrates its role in this paper, as a 
lower dimensional analogue of the Skyrme model. In fact Skyrme himself [12] proposed the sine- 
Gordon model as a toy model for a nonlinear meson field theory with the baryon interpretation 
of solitons ie the forerunner of the Skyrme model. Let us consider the Skyrme lagrangian 
(2.1) (including the pion mass term) not in (3+l)-dimensions with an SU(2)-valued field, but in 
(l-l-l)-dimensions with a U(l)-valued field. We can then write 

U = e''^ (2.6) 

where is a real-valued field, upon which (2.1) becomes (after dividing through by the constant 
4 to obtain the standard normalization) 

1 

^sg = ^d^,<f>d^'ct> + —icos<t>- 1) (2.7) 

where now x^^ = {t,x), fi = 0,l. This is of course the sine-Gordon lagrangian. Note that in three 
dimensional space it is not necessary to include the potential term (pion mass term) in order to 
have a finite size soliton solution, although the term quadratic in derivatives (the Skyrme term) 
is required. In contrast in one dimensional space the Skyrme term is not required (the Skyrme 
term contribution has vanished in (2.7)) but the potential term must be present. This can easily 
be seen as a consequence of the behaviour of the sigma model term (the first term) under a 
scale transformation of the space coordinates. In the intermediate case of two dimensional space 
both a Skyrme term and a potential term are required, and this produces the planar Skyrme-like 
model which we mentioned briefiy in the introduction. Since we require m 7^ in the sine-Gordon 
lagrangian we can choose length units to set m = 1 in (2.7) without loss of generality. 
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The sine-Gordon equation which foUows from (2.7) is 



d^(t)d'^(t) + sm(t) = 0. (2.8) 

The soliton solutions of the sine-Gordon equation (which are known as kinks) have finite energy 
so that the field ^{x) at fixed time must tend to an integer multiple of 27r as a; — > ±00. As with 
the Skyrme model there is therefore an integer- valued conserved topological charge Bgg, given by 

+ OC 

Bsg = ^ J d^dx. (2.9) 

— 00 

Again this counts the number of solitons (kinks) in a given configuration and we regard it as 
the analogue of the baryon number B of the Skyrme model. There is also a lower bound on the 
energy of a configuration in terms of the soliton number ie. E > 2'K\Bsg\. 

The static one-kink solution is given by 

(/> = 4arctan(e^-'') (2.10) 

where a is the position of the kink. This solution saturates the energy bound so that it has energy 
E = 2. Note that this is an important difference between the sine-Gordon model and the Skyrme 
model since it implies that there can be no soliton bound states in the sine-Gordon theory 
(although there are soliton- antisoliton bound states known as breathers). Another difference 
between the two models, which will be important in the next section, is that there is no internal 
phase or orientation for a kink, in contrast to the already mentioned orientation of the Skyrmion. 

There are no static two-kink solutions to the sine-Gordon equation, but the integrability of 
the theory allows the explicit construction of the exact two-kink solution where both kinks have 
arbitrary velocity (but at least one non-zero). The solution 

/■usinh(7a;) > _ ^ 

6 = 4arctan( , / : ) (2.11) 

^ ^ cosh(7ut) ' ^ ' 

where 7 = (1 — m^)~5,0<u<1, was first found by Perring and Skyrme [13] and describes two 
kinks which approach along the a;-axis with velocity u, scatter elastically at t = 0, and emerge 
from the interaction with a phase shift given by 

^=^i^^. (2.12) 

7 

This 2-kink solution has energy E = Aj. 

In the following section we briefly review the principal idea of the Atiyah-Manton approach 
to a finite dimensional truncation of the Skyrme theory, and then we construct in detail an 
analogous finite dimensional truncation of the sine-Gordon theory. 
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3. Kinks from (DP^ Instantons 



Let be the su(2)-valued Yang-Mills gauge potential for an instanton in four-dimensional 
euclidean space with coordinates = (x, 0:4), = 1,2,3,4. Note that the euclidean time X4 in 
the instanton theory is quite distinct from (and not to be confused with) the "physical time" t 
of the soliton theories introduced in the previous section. The main idea of the Atiyah-Manton 
scheme [1] is to construct the holonomy 



+00 



Z7(x)='Pexp J -A4{x,X4)dx4 (3.1) 

—00 

where V denotes path ordering. 

The fundamental topological result is that if wc consider this holonomy as a Skyrme field 
then it has baryon number B = n, where n is the topological charge (ie instanton number) of the 
gauge potential A^. The moduli space of charge n SU(2) Yang-Mills instantons has dimension 
8n — 3, and this generates a manifold of Skyrme fields with dimension 8n — 1 (after including 
a choice of gauge at infinity). It turns out that this manifold appears to be a good truncation 
of the Skyrme model. Here we will only mention the example for n = 1, where a 1-instanton 
positioned at the origin with scale A generates a Skyrmion of the hedgehog form (2.5) with a 
profile function 

This Skyrme field has minimum energy when A = 0.69 at which it takes the value Ei = 1.24x 127r^, 
which is within 1% of the numerically determined value Ei. 

There is a lower dimensional analogue of this procedure which was briefly introduced by the 
present author in [14]. Here we describe this result in more detail and use it to construct finite 
dimensional manifolds which we use in section 5 to study soliton dynamics. 

We use the gauge field formulation of the (DIP^ cr- model in two-dimensional euclidean space. 
It is defined in terms of a 2-component column vector Z, which is a function of the euclidean 
spacetime coordinates x^ = (x^,a;^), and is constrained to satisfy the condition 

Z^Z = 1. (3.3) 

The action density has a U(l) gauge symmetry and is given by 

£ = Tr{D^Z)'i{D''Z) (3.4) 
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where Tr denotes trace and are the covariant derivatives 

D^ = di,- (3.5) 
with the composite gauge fields being purely imaginary and defined by 

= Z^di^Z. (3.6) 
The equation of motion resulting from (3.4) is 

[a^W,P]=0 (3.7) 
where ]P is the one-dimensional hermitian projector 

P = ZZl (3.8) 
Instantons are finite action solutions to (3.7) and are most easily written using the parametrization 



^ ^1 + \W\^ [wj ^^'^^ 
where they are given by a rational function of x+ = xi + ix2- 

The target manifold is (DP^, which is isometric to 5^, so that due to the homotopy relation 

n2(5') = Z (3.10) 

each finite action field configuration has an associated integer winding number. This winding 
number is the topological charge (instanton number) of the configuration and is given by 

where Im denotes the imaginary part. For instanton fields the degree of the rational function 
W{x+) is equal to the instanton number n. 

Following (3.1) we construct the holonomy 

U{xi) = {-ire^p{J A2{xi,X2)dx2) (3.12) 

— oo 

which is a U(l)-valucd field. Note that there is no path ordering required here since we are 
dealing with an abclian gauge theory. The prefactor in (3.12) may be intcrpretated (see below) 
as the extra holonomy along a large semicircle at infinity in the {xi,X2) plane which connects the 
points X2 = +00 and X2 = — oo and arises from the fact that the instanton may be in a singular 
gauge at infinity. In fact the holonomy (3.1) also requires a prefactor of -1 if the instanton is in 
a singular gauge at infinity [1]. 
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We now make use of our interpretation of the sine-Gordon model as a U(l) Skyrme model 
through the relation (2.6) ie from the holonomy (3.12) we construct the sine-Gordon field through 
the definition 

The fundamental topological result is that the sine-Gordon field defined in this way has a kink 
number equal to the topological charge of the instanton gauge potential ie Bgg = n. This can 
be proved using the following steps. Combining (3.12) with (3.13) we have that 



-l-oo 



(x) = n(mod2)7r — i / A2dx2 



(3.14) 



Let L denote the line which joins the point {xi, — oo) to the point (a;i, -|-oo), and C be the large 
semicircle at infinity which joins these points with opposite sense of direction. Then we can write 
(3.14) as 



(j){x) = n(mod2)7r + i A^dx^ —i / A^dx^^ 



(3.15) 



The second term may be evaluated to give — titt which cancels exactly the first term (since <p 
takes values in we can freely add integer multiples of 2n). The remaining term is an integral 
over a closed contour so we can use Stokes' theorem to write it as an area integral. In terms of 
W the result is 

(l>{x) = J dxi j dx2 ^^^|^|2-)2 • (3-16) 

— oo — oo 

Comparing with equation (3.11) we see that (f){x) is exactly 27r times the instanton topological 
charge contained in the region of the {xi,X2) plane given by xi < x. Hence (j){—oo) = and 
4>{oo) = 27m, and we have the final result 

+00 

Bsg = ^ j d^<f> dx = ^{(l){oo) - (l>{-oo)) = n (3.17) 



The moduli space of charge n CP^ instantons has dimension An — 1. Since the holonomy 
(3.12) is invariant under translations of X2 the above procedure produces kink fields which depend 
on a manifold, which we shall call M.^, that has dimension 4n — 2. 
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Let us consider the case n = 1. Then we may take W to be given by 



W = X{x+-a) (3.18) 

where A is real and positive and is the scale of the instanton, and a may be taken to be real due 
to the above mentioned invariance. This generates the 1-kink field 

= cl>nx, A, a) = n{l + —^=£==). (3.19) 

^/l + A^(a; — a)^ 

We therefore have that Ad^^ = R"*" x R, with the interpretation of the coordinates as a scale and 
position for the kink. The energy of (3.19) is obviously independent of a, and is minimized when 
A = Acp = 0.695 at which it takes the value E = Mcp = 2 x 1.010. This is the mass of the (CP^ 
generated approximate kink, which is within 1% of the energy of the exact 1-kink solution (2.10). 
Let us note that not only is the procedure wo have described very similar to the Atiyah-Manton 
procedure but the results in the 1-soliton sector have a remarkable similarity also. Compare the 
form of the approximate kink (3.19) with the approximate profile function (3.2), the values of 
the instanton scales at which the approximate soliton energies are minimized, and the percentage 
error by which the energy of each approximate soliton exceeds that of the exact solution. 

Now consider the case n = 2. In order to make the interpretation of each coordinate as clear 
as possible we write W in the form 

= r + , ' . , (3.20) 

(x+ - aij {x+ - a2 + ic) 

where Ai,A2 are the positive scales of each instanton and ai,a2,c are real and give the position 
on the xi axis of each instanton and their separation in X2 space respectively, and 6 gives the 
relative internal phase between the two instantons. 

We have that M^" = x x (R"*")^. However, as we shall demonstrate shortly, the 
manifold A^2^ is not a good truncation of the 2-kink sector, unless we restrict to a certain 
submanifold. Before we begin this discussion let us make a few simplifications by considering 
only two-kink configurations in which the kinks are placed symmetrically about the origin and 
have momenta of equal magnitude but opposite sign ie we arc considering 2-kink solutions that are 
given by the exact solution (2.11). This corresponds to setting Ai = A2 = A and ai = —02 = a; we 
also choose to set c = for simplicity. Our first argument is no more than an intuitive observation 
as to why the value of 6 should be taken to be zero. We then make an explicit computation to 
show that the manifold obtained by allowing 9 to take arbitrary values results in an incorrect 
approximation to the sine-Gordon model, since it predicts a spurious 2-kink bound state. 
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As noted earlier an important difference between tfie Sfcyrme theory and the sine-Gordon 
theory is that the sohtons of the latter have no internal degrees of freedom ie no orientation or 
phase. We therefore conjecture that to approximate sine-Gordon solitons one should, in contrast 
to the Skyrme case, consider only instantons which are all in phase. For n = 2 this means setting 
^ = 0. We have now restricted to a two dimensional submanifold of M!^ which we denote by 
M!^ ■ In fact A^2^ = R"*" x M"*", where the coordinates are a and A. Note that a is no longer 
allowed to be equal to zero since the two-instanton solution would then degenerate into a one- 
instanton solution. Explicitly from M!^^ we obtain (after calculating the integral in (3.12) by 
contour integration) the 2-kink field 



1 

^a;2 + o2 + 2A-2 + 2^2(3.2 + ;^-2) + ^ + ^3.2 + ^2 + 2A-2 - 2Va2(a;2 + A-^) + A-^ 

(3.21) 

If we compute numerically the energy of this configuration as a function of a and A we find that 
it is minimized in the limit o ^ 00, in which (3.21) becomes asymptotically two approximate 
kinks of the form (3.19) with positions a and —a. The energy in this limit is therefore minimized 
when A = Acp at which it takes the value E = IM^p. M.^ therefore produces the correct 
qualitative sine-Gordon behaviour of no bound states in the two-soliton sector. Note that this 
feature is not automatic for our approximate 2-soliton sector since the approximate 1-soliton does 
not attain the topological lower bound ie Mcp > 2. We will now demonstrate that it is precisely 
this qualitative "no bound states" requirement which is violated by Al2^ as a finite dimensional 
truncation. In the same way that we defined the submanifold we define the similar two 
dimensional submanifold M.2' but where we set ^ = tt rather than the previous case of 6 = 0. 
We are now considering the case in which the two instantons are exactly out of phase. Note that 
■^2^ = IR X since the position coordinate a is now allowed to take the value zero, since this 
2-instanton solution no longer degenerates into a 1-instanton solution in the a ^ limit now 
that the instantons are not identical (they have opposite orientation). Explicitly from A^2^ we 
obtain the two kink field (where for notational convenience we define C = A2a;2 — 1) 
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02 A,a) 



nx 



VC-Xa 



< ^ y/X{x^+a'^)+2aVC ^/\{x'2+a^)-2aV< 



ifC>0 



■Kxav2\ 



X 



(3.22) 



{C-\/ v^A2(^ + a2)2 + 4a2c2 + A(x2 + a2) 
^ +Aa-^V'A2p + a2)2 + 4a2C2 - A(x2 + a2)} if C < 0. 



Of course in the limit a ^ oo this solution again becomes two approximate 1-kinks as before, 
with the energy minimized at A = Acp and taking the value 2Mcp. However, the crucial point 
is that the energy of (3.22) is not minimized by this limit of a. Numerically calculating the 
energy of (3.22) we find it is minimized at a = 3.16 and A = 0.55, at which it has the value 
E = 4 X 1.003 < 2Mcp. The exact values of these constants are unimportant, the point being 
that the energy is minimized at a finite value of a. This means that M.2' does not provide a 
qualitatively good truncation of the 2-soliton sector since it describes a theory with a 2-kink 
bound state. This observation is a warning that care must be taken when applying the instanton 
method to the study of solitons. If we had been studying a non-integrable soliton theory, with no 
exact results to guide us, then the instanton method may have led us to the incorrect conclusion 
that a 2-kink bound state existed in the full theory. It would appear that the only way to 
attempt to avoid making such errors in non-integrable theories is the kind of intuitive argument 
given above, based upon relating the instanton orientation to the soliton orientation (or lack of 
it). 

It is in principal possible to calculate explicitly the 2-kink field corresponding to the manifold 
A42', since using contour integration methods it depends on the zeros of a quartic polynomial, for 
which a closed form expression exists. However, we have already demonstrated that a submanifold 
of A^2^, namely does not provide a good truncation and this implies that neither does ^^2^. 
This supports our conjecture that the correct manifold to use is A^2^ obtained from instantons 
with the same orientation. In the rest of this paper we will restrict to this manifold when 
discussing the 2-soliton sector generated from (DP^ instantons. 
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4. Kinks from Yang- Mills Instantons 



As briefly described earlier the Atiyah-Manton procedure involves taking SU(2) Yang-Mills 
instantons in R"* and computing a holonomy. In the previous section we introduced an analogue 
of this procedure which again involves computing holonomics but this time of (DIP^ instantons in 
R^. There is yet another analogue, which we shall now discuss, where we begin with the same 
instantons used in the Atiyah-Manton scheme. 

It is obvious that to obtain sinc-Gordon kinks from SU(2) Yang-Mills instantons in R'* one 
has to do something other than simply compute a holonomy. The problem is how to construct 
a real- valued field in R from such an instanton. The solution is suggested by equation (3.16) of 
the previous section. There we saw that computing the holonomy (3.12) was equivalent, in that 
case, to computing the integral of the topological charge density over a certain region of space. 
For other instanton systems there will be no such equivalence but it is clear that the second 
interpretation can be taken as a definition to define the procedure. In reference [15] the authors 
considered such a procedure by defining 



-|-oo -|-oo 



(t){x)=2TT j dxi j dx2... j dxd qd{xi,...,Xci) (4.1) 

— oo — oo — oo 

where qd denotes the topological charge density of any instanton theory in d-dimensional euclidean 
space. Clearly by construction if wc consider an n-instanton solution to the gauge theory in 
question then (p{oo) = 2Tm and 0(0) = so that Bgg = n. This procedure becomes the holonomy 
procedure of the previous section if we choose d — 2 and the instanton theory to be the (CIP'^ a- 
model. The analysis in [15] was concerned only with the 1-soliton sector where it was shown that, 
for the theories they considered, the best approximate 1-kink was obtained from four-dimensional 
SU(2) Yang-Mills theory. In this section we will therefore concentrate on this particular instanton 
theory. 

Let be the su(2)-valued Yang-Mills gauge potential for an instanton in four-dimensional 
euclidean spcice with, coordinates x^^,, /x — 1,2,3,4. The gauge field is F^^, = d^A„ - d^A^ + 
[^^, Ay]. In terms of these quantities the formula (4.1) is 



-t-oo 



<^(a;) = Y^ / '^'^i / dx2dxzdxi e^,^pTT{Fi"' F'^^) (4.2) 

— oo — oo 

where is the totally antisymmetric symbol on four indices. 

The moduli space of charge n SU(2) Yang-Mills instantons has dimension 8n — 3, but since 
the above integral is invariant under shifts of Xj for j = 2, 3, 4, this procedure produces kink 
fields which depend on a manifold, which we call M.'^, which has dimension 8n — 6. 
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One could now simply substitute an n-instanton solution into (4.2) and attempt to perform 
the integration, as done in [15] for the n = 1 case. However, let us first comment that, as 
mentioned above, (4.2) is to be thought of as a similar expression to (3.16). It would be interesting, 
and certainly extremely useful from a computational point of view, if another form of the relation 
(4.2) could be found which was closer in spirit to the form (3.12); in the sense that an integral 
over xi did not have to be performed. In the following we restrict to instantons obtained from 
a superpotential [16, 17] (ie 't Hooft or the more general JNR instantons) and show that such a 
relation exists. Note that for the cases n = 1, 2 such a restriction results in no loss of generality 
since all instantons arc of this form. It is also possible to produce a similar result in the general 
instanton case, using the ADHM [18] construction and an expression for the topological charge 
density as a double laplacian [19], but we shall not describe this here. 

The instantons we consider are therefore given by 

A^='-a^,d''\ogp (4.3) 

where cr^^ — {^i^n/a + S^^aSii, — S/^f^S^a)'^"' , a = 1, 2, 3 are the Pauli matrices and /^(x^) is the 
real-valued superpotential satisfying the wave equation in euclidean four-space Op = 0. Then 
(4.2) becomes 

X +CXD 

<p{x) = J dxi J dx2dx3dx4\3\3 logp (4.4) 

— oo — oo 

where, as in [17], this integral must be interpreted as excluding infinitesimal regions of integration 
containing the singularities of p; this is equivalent to defining □ □log(a;^ — X^)'^ to be zero for 
any constant in ]R^. By converting total derivative terms to boundary terms we write (4.4) 
in its final form 

-t-oo 



x) = n(mod2)7r — / dx2dx3dx4 dflog p 



(4.5) 



where we have succeeded in removing the integral over xi . 

Take the case n = 1, where we may write (using the above mentioned shift invariances) 

P = '^ + l ^ \ ,2 - (4-6) 

(x^ - adi^y 

Here A is the positive scale of the instanton and a is real, giving the position of the instanton 
on the xi axis. A4Y^ = R"*" x ]R and corresponds to the explicit 1-kink field (obtained from 



14 



integrating (4.5)) 



rn^, A, a) = TT + ^dl{{x - aWX-^ + {x-ar } 

A(x-a) lX{x-a) (4-7) 

Vl + A2(a;-a)2 (1 + A2(ar - a)2)i ^ 

The final expression is the form in wliieh it was obtained in [15] by calculating explicitly the 
integral (4.2). The energy of (4.7), independent of the kink position a, is minimized when 
A = Xym = 0.419 at which it takes the value E = Mym = 2 x 1.0002. This is the mass of 
the Yang- Mills generated approximate kink, which exceeds that of the exact solution (2.10) by 
around 0.02%. It is therefore a remarkably accurate approximation. 

We now consider the n = 2 sector, and following the remarks of the previous section, consider 
instantons which have no relative orientation. We have not performed an explicit calculation to 
see whether the "no bound states" requirement is violated for Yang-Mills instantons with a 
relative orientation. Note that since {Mym — 2) is two orders of magnitude smaller than (Mcp — 2) 
the window for such a possible violation is greatly reduced. Accurate numerical work is required 
to settle this issiie. Wc also make the same symmetric kink simplifications that we made in the 
previous section, and consider the submanifold 7V(|™ = R"'" x J^,"*" corresponding to the instanton 
superpotential 

This gives the 2-kink field 

x{x^ + 3a2 + A-2 + ^(x2 - a2)2 + 2A-2(x2 + a2) ) ^ (4.9) 

= TVO^{ — I 

Va;2 + a^ + A-2 + y/I^^ + A-^ + ^ x^ + a2 + A-2 - V4a;2a2 + 

Note that upon performing the double differential in the above expression the result will be quite 
lengthy and would have been tedious to calculate from the expression (4.2) directly. 

The energy of (4.9) is minimized in the limit a — > oo, when the field becomes two approximate 
kinks of the form (4.7) with positions a and —a, so that the energy is minimized when A = Xym 
at which it is equal to 2Mym- 
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5. Dynamics on Ai'ji and 

In this section we compute the approximate sohton dynamics, in the one and two-soHton 
sectors, obtained from the finite dimensional truncations constructed in the previous two sections. 

We first consider Ad^^, with parameters a and A. As described in the introduction, time 
dependence is introduced by allowing these parameters to be functions of the "physical time" 
t ie a{t) and X{t). A lagrangian for this system is then obtained by restricting the field theory 
lagrangian to the configuration space described by A^^^ ie the lagrangian is obtained by substi- 
tuting the ansatz (3.19) into the lagrangian density (2.7) and integrating over x. Before we give 
the result it is convenient to introduce the normalized scale /i — X/Xcp, so that the minimum 
energy of a static configuration is obtained at the value /U = 1. The result for the lagrangian is 



/:f = ^(Ma^ + fccp4 (5.1) 



where dot denotes difFcrcntiation with respect to time t and wo have introduced the constant 
kcp = ~ 0.69. It is easy to see that the equations which follow from (5.1) are 



a/i + d/i = 

2kcpjj,iJ, — 2>kcpiJ? + — a^) — ji^ =0 



(5.2) 



and have a one parameter family of solutions (independent of kcp) given by 

a = ut, M = 7 (5.3) 

for — 1 < u < 1 constant, where 7 is the Lorentz factor {1 — u'^)~i . The energy of such a solution 
is "fMcp. Hence wc obtain a simple Lorentz boost of the static approximate kink. Note that 
to be able to obtain an exact Lorentz boost from the instanton scale is a feature unique to the 
approximation of solitons in one dimensional space. 

If we now consider Alf™, by performing the same procedure with the ansatz (4.7) wc obtain 

the same resiilts as above after a labelling replacement cp ym, where the only quantity we 

have yet to define is kym = w^— ~ 0.81. It is interesting to compare these two results with that 

obtained by taking the exact kink solution (2.10) as an ansatz and introducing a scale through 

the replacement x — > fix. We then get a similar result with the obvious notation Mgxact = 2 and 
2 

kexact = j2 ~ 0.82. This demonstrates that the coupling of translational and scale degrees of 
freedom is very similar for the approximate kinks as it is for the exact kink; compare kcp with 
kexact and particularly kym with kexact- 
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Let us now turn to the 2-soliton sector and consider M.^ . We are going to make one 
further simphfication (it turns out that we obtain almost identical results with or without this 
simplification). In the limit of large a (ie well-separated kinks) a 2-kink configuration becomes 
two 1-kink configurations so that the lagrangian becomes w ILf . The solution in this region 
will therefore be two free, Lorentz boosted, kinks with d = — u, A = 7Acp, corresponding to 
studying kink scattering with kink velocity < u < 1; compare the exact solution (2.11). Our 
simplification is to assume that the scale changes little during the evolution so that we can fix it 
at this value ie we set A = 0, and truncate to a system with only one degree of freedom. (This 
is in fact a good approximation. Wo have studied the two-dimensional system, which arises by 
allowing A(i), numerically using a Runge-Kutta evolution and find that the variation of A from 
the constant value A = ^\cp is so small as to have a negligible effect on our results). 

Substituting the ansatz </)2^(a;, 7Acp, a(i)) into (2.7) and integrating over x gives 



-2 



gcp{a)c? - Yap{a) (5.4) 



where the metric Qcp and the potential Vcp are given by 
9cp{a) = ^Jdx i^^{x,a,'yXcp)f 





oo 



(5.5) 



Vcp{a) = I dx -{-^{x,a,j\cp)) + -{I - cos(j)2^{x,a,jXcp)). 







Of course, in the Yang-Mills case wc have exactly the same formalism, again simply replacing 
cp by ym so that wc arc then using the fmiction (/)|™(a;, 7Aj,„i, a(t)) etc. In Fig. 1. wc plot the 
metric and potential functions ^cp, S'j/m) ^ym obtained by integrating the above expressions 
numerically. 
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Figl : The metric and potential functions gcp,Vcp (continuous curves) and gym,Vym (dashed 
curves). In each case the potential curves are the monotonic decreasing functions of a. 

To produce the plots of Fig 1. we have set 7 = 1. The range of a reflects the values which 
are relevant in the scattering processes we shall discuss below. It is easy to see that since is 
a symmetric function of a then Qcp — > as a ^ 0. Also it can be shown that the slope of at 
the origin behaves like so that the potential grows without bound as a — > 0. Similar remarks 
apply for 02"*- F^om (5.4) the equations of motion are 

with initial conditions a(t = 0) = ao, a{t = 0) = —u. 

The solution can be obtained implicitly by quadrature as 



where E = ^gcp{ao)u^ + Kp(ao) is the total energy. This solution is valid upto t = ti, where ti is 
the turning time ti = t{ai) and ai is the turning point Vcp{ai) = E. The position for t> ti can 
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be determined from the symmetric property of the motion a{t — ti) 
2-kink phase shift is given by 

Sep = 2(ao - uti) 



= a{ti —t). The approximate 
(5.8) 



in the hmit in which ao — > oo with ti determined as above. 

We now compare the results of our finite dimensional truncations with the exact results 
given in section 2. First note that to compare with the exact 2-kink solution (2.11) it is necessary 
to shift the time coordinate t in this solution to t — — 1 — oq), so that the kinks have 
positions ±ao at t = 0. In Fig. 2a. wc plot the fimction (f)2^{x,j\cp, a{t)) (continuous curve) and 
the exact solution (2.11) (dashed curve) for a velocity u = 0.5, at various times. Wc have taken 
flo = 10 and all required integrals arc performed numerically. Fig. 2b. is as Fig. 2a. but for the 
Yang-Mills function <p2"^{x,jXym,a(t)) with the time evolution of a determined by £2™- 
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Fig 2a : (continuous curve) and the exact solution (dashed curve), at increasing times. 
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Fig 2b : As Fig 2a. but for the approximation 
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Fig 2a. (1) shows the initial (t = 0) configuration, and it can be seen that it is very close to 
the initial exact solution. The kinks then collide, and Fig 2a. (2) shows the configurations at the 
time for which the separation a is a minimum. Note that since we are considering the scattering 
of identical kinks it is impossible to distinguish forward from backward scattering. Our instanton 
moduli coordinates make backward scattering the most natural interpretation (since a G ]R,^) 
but we interpret the results in the more usual forward scattering manner in order to compare 
with the exact results. We can see that even when the kinks have lost their individual identities 
the approximation is still quite good, although it is clear that there are two regions in which the 
error is now appreciable. Fig 2a. (3) corresponds to the time for which the position a has returned 
to its initial value. The error here is quite clear, with the approximation giving a phase shift in 
excess of the exact result. 

We now turn to Fig 2b. These figures correspond to those of Fig 2a. but now for the 
Yang-Mills generated approximation. Again two curves are graphed in each window, but the 
approximation is so close to the exact solution that it appears there is only a single graph! This 
approximation is clearly extremely accurate and much the superior of the two schemes. This is 
to be expected since the one-kink field ^2"* is a much better approximation than the one-kink 
field <f)l''. 

In Fig. 3. we plot the exact phase shift S (continuous curve), the (DP^ generated approxi- 
mation Sep (stars) and the Yang-Mills generated approximation Sym (circles) for velocities in the 
range 0.1 < u < 0.9. This figure displays the main result of this paper, showing how accurately 
the instanton truncations approximate the infinite dimensional field theory dynamics. Although 
the cp approximation gives qualitatively correct results the errors are quite large. This may seem 
surprising when one considers the accuracy of the one- kink field recall its energy is only 1% 
larger than that of the exact one- kink. However, an important point to bear in mind is that the 
approximation decays power-like, as opposed to the exponential decay of the exact one-kink. This 
would therefore suggest that we would be over optimistic to hope to obtain results of too great 
an accuracy. This is where the ym approximation yields surprising results. It too generates kinks 
with a power-like tail behaviour, but from Fig. 3. we see that the accuracy is remarkable, for 
the whole range of velocities. So it would seem that the exponential localization of sine-Gordon 
solitons is not a major factor in influencing their interaction. This is an interesting point which 
requires further investigation and is of relevance to higher dimensional systems where solitons 
are often localized in a power-like way eg lumps in planar sigma models. 



22 




23 



6. Conclusion 



This paper has dealt with some simpler analogues of the Atiyah-Manton construction of 
Skyrmions from instantons. We have shown how their procedure may be adapted to study other 
topological solitons, and furthermore have demonstrated the validity of the scheme in a case 
where exact results arc available for comparison. The discovery that constraints may have to be 
imposed upon the instanton orientations, in order to correctly reproduce the qualitative bound 
states behaviour of the soliton theory, is unexpected. It is clearly an important issue to consider 
when applying the scheme to non-integrable soliton theories. For the particular case of the 
Skyrme model the intuitive argument described earlier suggests that no constraint is required, 
since the orientation of the Yang-Mills instanton produces the required isospin of the Skyrmion. 

The results in this paper are encouraging for the use of the instanton method in approxi- 
mating the dynamics of Skyrmions. As mentioned earlier, there are as yet no complete results 
on the instanton truncation of Skyrmion dynamics, mainly due to the heavy computations this 
involves, and some technical difficulties in computing the non-abelian holonomy. Using numeri- 
cal methods [20] there are some results regarding the potential function on the finite dimensional 
two-Skyrmion manifold, and perhaps a full numerical scheme could also be used to construct the 
metric function and solve the resulting ordinary differential equations which result. 

In this paper we have dealt with only the classical dynamics of solitons. It is also possible 
to study quantum scattering by a quantization of the finite dimensional system obtained from 
instanton moduli. A finite dimensional truncation of the sine-Gordon model (obtained from 
patched kinks) was quantized in [21] and found to compare well with the known exact results 
in the weak coupling limit. A similar analysis carries over to the finite dimensional manifolds 
constructed in this paper and results of a similar quality can be obtained. 

It appears that the soliton from instanton construction is a useful approximation technique, 
the wide applicability of which is only just beginning to emerge. It is already known that this 
scheme can be applied to other kink systems, such as the double sine-Gordon equation [15], to 
produce accurate approximate static kinks, and it produces excellent results when applied to (fP'^ 
field theory [22] . No doubt the method will prove to be a useful tool as other areas of application 
come to light. 

Finally, let us make a few remarks about Skyrmions which arise from the work in this paper. 
One point to note is that our earlier comments about the power-like decay properties of the 
approximate kinks, as compared to the exponential decay of the exact kink, also have relevance 
in the Skyrme context. The instanton generated profile function (3.2) decays like ^ for large r, 
as does the (numerical) exact solution. However, if a pion mass term is added then the exact 
solution now has an exponential decay. The results of this paper would seem to suggest that 
this is not a serious defect if one wished to use the Atiyah-Manton scheme to study Skyrmion 
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dynamics even in the massive pion case. A further observation is based on the similarity of the 
Skyrmion profile functions (3.2) and the kink function (3.19). In fact, if one considers a (D]P^ 
anti-instanton (rather than an instanton) positioned at the origin and makes the replacement 
a; — > r then one obtains exactly the profile function (3.2). This identification prompted the use of 
a sine-Gordon anti-kink field as a trial function for the Skyrme profile function [23] and slightly 
reduced the energy as compared to (3.2). Now note that we can use a similar identification with a 
Yang-Mills generated anti-kink, to propose a modification of the Atiyah-Manton profile function 
(3.2) to 

/(r) = 7r(l - , ^^^). (6.1) 

' ^ Vl + A2r2 (1 + A2r2)i^ ^ 

This function is the same as (3.2) but with an additional term added which docs not affect the 
boundary conditions. As described earlier (3.2) gives a Skyrme field with energy about 1% above 
the (numerical) exact solution, whereas (6.1) reduces this excess to around i%. If one replaces the 
^ coeffecient of the final term in (6.1) with an arbitrary constant c, then the energy is minimized 
for c « 0.4 at which the excess energy is less than By introducing independent scales for 

the two terms in (6.1) one could perhaps reduce the energy still further. 

Acknowledgements: I thank the EPSRC for a research fellowship. 



25 



REFERENCES 



1. M. F. Atiyah and N. S. Manton, Phys. Lett. 222B(1989), 438. 

M. F. Atiyah and N. S. Manton, Comm. Math. Pbys. 153(1993), 391. 
N. S. Manton, ^^Skyrme fields and instantond^ 

in Geometry of Low Dimensional Manifolds: 1 (LMS Lecture Notes 150), 
eds. S.K.Donaldson and C.B.Thomas: Cambridge University Press 1990. 

2. N.S. Manton, Phys. Rev. Lett. 60B(1988), 1916. 

3. M. F. Atiyah and N. J. Hitchin, " The Geometry and Dynamics of Magnetic 
Monopoles " (Princeton University Press, Princeton, 1988). 

4. T.M.Samols, Comm. Math. Phys. 145(1992), 149. 
P. J.Ruback, Nucl. Pbys. 296B(1988), 669. 

L A. B. Strachan, J. Math. Pbys. 33(1992), 102. 

5. R.S.Ward, Phys. Lett. 158B(1985), 424. 
R.A.Leese, Nucl. Phys. 344B(1990), 33. 

6. G.W. Gibbons and N.S. Manton, Nucl. Phys. 274B(1986), 183. 
B. J. Schroers, Nucl. Phys. 367B(1991), 177. 

7. P. M. Sutcliffe, JVonJinearitj 4(1991), 1109. 

8. P. A. Shah, "Vortex scattering at near critical coupling' , 
Cambridge preprint DAMTP 94-8. 

9. T. H. R. Skyrme, JVuci. Piiys. 31(1962), 556. 

10. E.Witten, Nucl. Pbys. 160B(1979), 57. 

11. G.S.Adkins, C.R.Nappi and E.Witten, JVucJ. Piiys. 228B(1983), 552. 

12. T. H. R. Skyrme, Proc. R. Soc. 262A(1961), 237. 

13. J. K. Perring and T. H. R. Skyrme, Nucl. Phys. 31(1962), 550. 

14. P. M. Sutcliffe, Pbys. Lett. 283B(1992), 85. 

15. G. N. Stratopoulos and W. J. Zakrzewski, Z. Phys. C 59(1993), 307. 

16. G. 't Hooft, unpublished. 

E. F. Corrigan and D. B. Fairlie, Phys. Lett. 67(1977), 69. 

F. Wilczck, "Quark confinement and field theory" 

ed. D. Stump and D. Weingarten (John Wiley, New York, 1977). 

17. R. Jackiw, C.Nohl and C.Rebbi, Phys. Rev. 15D(1977), 1642. 

26 



18. M.F.Atiyah, V. G. Drinfield, N.J. Hitchin and Yu. I. Manin, Phys. Lett. 65A(1978), 
185. 

19. E. F. Corrigan, D. I. Olive, D. B. Fairlic and J. Nuyts, Nucl. Phys. 140B(1978), 31. 

20. A. Hosaka, S. M. Griffics, M. Oka and R. D. Amado, Phys. Lett. 251B(1990), 1. 
A. Hosaka, M. Oka and R. D. Amado, Nucl. Phys. 530A(1991), 507. 

21. P. M. Sutcliffe, iVud. Piiys. 393B(1993), 211. 

22. P. M. Sutcliffe, Plijs. Lett. 328B(1994), 84. 

23. P. M. Sutcliffe, Plijs. Lett. 292B(1992), 104. 



27 



